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•The effect of noise on a rotational mode of a pendulum excited kinematically in vertical direction 
as been analyzed. We have shown that for a weak noise transitions from oscillations to rotations 
* ^a. nd vice versa are possible. For a moderate noise level dynamics of the system is governed by a 
recombination of the excitation amplitude and stochastic component. Consequently for stronger noise 
l_Jhe rotational solution as an independent sychonized mode has vanished. 
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"^I" Introduction 

o ■ 

forced pendulum is one of the simplest nonlinear system which shows typical 
^Ochaotic behaviour and this has been discussed extensively in many research 

particles and in several textbooks, for example see [1-11]- It is worth noting the 

Cpendula have been studied not only to understand the fundamental nonlinear 
>-^behaviour but also with a view to solve some practical problems of wave energy 

Hextraction from sea waves [12], gravitational gradient pendulum [13,14] and 
.^hip dynamics [15], to name a few. Some other nonlinear phenomena like the 
^Josephson junction [16,17] and can also be easily associated with the analysed 

S-jproblem. 

■ - - In this note we will examine stability of rotational motion of parametric 
pendulum [6,11] in noisy conditions [18]. For a pendulum under the harmonic 
kinematic excitation (Fig. the equation of motion can be written as follows 



(1) 



2 



Phase lock and rotational motion of a parametric pendulum 



where m is point mass I is a length of pendulum, g is a gravitational constant. 
Introducing the natural frequency of free oscillations, ujq = yj g/l we define 
dimensionless time r = oj^t and renormalized frequency £1 £1 = fi'/ w o- 
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Figure 1. (a) A parametric pendulum with harmonic vertical excitation of a suspension point 
x = Acosf2't = Acosf2r, (b) Rotational number < v > /Q versus 7 for a = 0.1 (v = tp). Vertical 
lines correspond to 7 = 0.25 and 0.50 (O = 2), respectively. Calculations have been done starting 
from the smallest 7 fa and increasing it up to 7 = 0.9. For each new 7 the initial conditions <j>o 
and vq were defined by the final values of <f> and v for previous 7. Numbers: '1', '2', '3' and '4' are 
related to synchronized motions: clockwise rotation, oscillation, anti-clockwise rotation and chaotic 

motion, respectively. 




Figure 2. Phase portraits and Poincare maps of the pendulum (Eq. 2) computed for u) = 2 
7 = 0.25 and 7 = 0.50. In Fig. 12a numbers '1', '2' and '3' correspond to solutions marked in Fig. [Tb 

(H = 2). 
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Figure 3. Rotation number < v > /fi (Q = 2) versus noise level a for a = 0.1 and four values of 
Yj (a) 7 = 0.4, (b) 7 = 0.6, (c) 7 = 0.8, (d) 7 = 1.0. 



Note the time dependent term ^j- cos $7t sin <fi is directly related to the iner- 
tial force due to the kinematic harmonic excitation with 0' frequency (Fig.^i). 



x(t) = A cos Q't = A cos fir. 



(2) 



To simplify the notation and further analysis we also introduced two new 
parameters 7 = A/ 1 and @ = k/(ml 2 ) while time derivative d/dr = ' 



+ (54> + (1 + 7ft 2 cos fir) sin = 0. 



(3) 



Deterministic system regular oscillatory and rotational motions as well as 
chaotic oscillations can be easily classified by observing changes in so called 
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Figure 4. Phase portraits and Poincare maps for a = 0.1, a = 0.05 and four values of 7; (a) 
7 = 0.4, (b) 7 = 0.6, (c) 7 = 0.8, (d) 7 = 1.0. Note different scales in «-axis. 

rotational number < v > /Q denned in [1, 19] as 

1 f T2 ■ 

< v > /n = lim — / 4> dr, (4) 

where T\ was the time interval necessary for the system to reach a steady state 
(here Ti = 2.3 x 10 3 ) while T2 was chosen to be large enoug h (t 2 = 4.3 xlO 3 ). In 
Fig. ^> we show that this quantity versus 7 and one can clearly see the regions 
of synchronized motions and their transitions into chaotic oscillations. Note 
for the synchronized motion the averaged rotational velocity is kept constant 
against 7 parameter up to 7 ~ 0.47 where we observe different behaviour (Fig. 
^>). Namely < v > /Q = ±1 for rotations and for an oscillatory motion. The 
rational number of < v > /O, exhibits the phase lock phenomenon [1] charac- 
teristic for nonlinear systems. The simulations of deterministic equation (Eq. 
show that in the relatively large region of system parameters the average 
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Figure 5. Time series showing the influence of the excitation strength, 7, computed for a = 0.1, 
a = 0.05, = 2; (a)7 = 0.4, (b)<y = 0.6, (0)7 = 0.8, (d)7 = 1.0. Note different scales in u-asds. 

velocity is constant. Consequently pendulum motion can be decoupled into 
two modes: rotations with the average rotational velocity and superimposed 
oscillations. To illustrate changes in corresponding attractors we plotted (in 
Fig. and b) phase portraits for 7 = 0.25 (and various initial conditions) 
and 7 = 0.5 representing different types of motions. 

In the next section we will examine the effect of noise in the examined 
system (Eq. |3J). In particular, we will discuss the simulation results of the 
rotational number in chaotic and noisy conditions, focusing on destabilization 
of rotational motions. 
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Figure 6. Rotational number < v > /Q (Q = 2) versus 7 for a = 0.1 and two different noise levels; 

(a) cr = 0.01, (b) a = 0.1. 



2 Phase lock in presence of noise 



To examine the system further we have transformed the second order differ- 
ential equation into two equations of the first order 

= v (5) 
ii = —f34> — (1 + 7O 2 cos Or) sin <f>. 

We assumed that in the system effected by noise through a kinematic exci- 
tation is not perfectly periodic resulting in variations of forcing frequency Q' 
or its related period 2tt/Q' in time. Having this in mind we decided to mimic 
this effect by a bounded noise concept [20] This type of noise is introduced by 
a random phase to kinematic forcing term cos (Or + there = ^(r) is a 
Wiener process. To perform numerical calculations we have to discretise our 
system (Eq. |3J)- First we write the Euler-Maruyama scheme of integration [21] 

v n+1 = (-(3v n - (1 + 7O 2 cos 6„) sin 4> n ) At 

(pn+l = "nAr (6) 

e n+ i = &n + OAr + T(r, cj)v / A7, 

where T(t) is a random number with a normal distribution and the standard 
deviation a referred here as the level of noise. To improve precision of our 
numerical calculations we have used the Runge-Kutta-Maruyama algorithm 

[21] which treats more carefully the deterministic part of Eq. H3 Figures 
l3K~d show the rotational numbers < v > /Vt versus a for different 7. Note 
that in all cases rotational motion is destructed by noise, and large enough a 
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the resulting cases have variable average velocity. In Fig. the deterministic 
state is originally chaotic for deterministic case and hence does not change 
in noisy conditions. Namely the result correspond to the chaotic oscillations 
with a noise component. For better clarity we have also plotted the phase 
portraits and Poincare maps for one the chosen value of a = 0.5 and the same 
values of 7 (7 = 0.4, 0.6, 0.8, and 1.0 for Fig. 0^,-d respectively). Note, first 
three figures (Fig. 0^,-c) show the system with large amplitude variations but 
not synchronized rotations; this appears to be chaotic oscillations without any 
clear synchronized rotation < v > /O ~ 0. On the other hand Fig. 0Ji is rep- 
resenting clock-wise rotations with the fixed average < v > /O = — 1. The 
corresponding time series are presented in Figs.EH-d. Phase lock phenomenon 
is present in Fig. EJl where we observe anti-clockwise rotation (persistent neg- 
ative velocity sign). Surprisingly the succeeding Poincare points are smeared 
by random excitations leaving average velocity < v > /O = — 1 unchanged. In 
Fig. one can see intermittent rotation with a negative sign. Note that in 
this case for a = the system behaves regularly (Fig. |2K) with clear clockwise 
or anti-clockwise rotations ( < v > /Q. = ±1). Consequently one can identify 
a path to chaotic oscillations induced by noise. Finally, in Figs. and b we 
show < v > /Q, versus 7 for two different levels of noise; a = 0.01 for Fig. 
EH and a = 0.1 for Fig. Et>, respectively. For a = 0.01 (Fig. EH) we observe a 
similar sequence of transitions from oscillatory to rotational motions as in the 
deterministic case (Fig. [Dp), while the noise level, a = 0.1, as shown in Fig. 
|U>, is strong enough to destroy the synchronized rotations for any 7. 



3 Summary and Conclusions 

We have performed numerical analysis of a parametrically excited pendulum 
focussing on the rotational motion and the influence of noise in the chaotic 
regimes. The rotational motion has been computed for the pure deterministic 
system and for the case with weak bounded noise introduced through a random 
phase in the excitation term. We found out, that in a large interval of increasing 
noise level the average rotational velocity was stable; the noise component 
created oscillations around the rotational average velocity. However, for a large 
enough noise level a = 0.1 (Fig. EH) the rotational motion as an independent 
sychonized mode has vanished. 

Looking for a rotational motion of pendulum we observed an intermittent 
transition to chaotic motion induced by noise. Such a conclusion was also given 
by Blackburn in his recent paper [18] where the effect of potential well escape 
was analyzed assuming a random external force. In our case motivated by a 
possibility of sea wave energy extraction by a parametric pendulum device, we 
considered stochastic distribution of the wave phase in succeeding periods. 
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